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Abstract 

Electron instabilities in the Hubbard model with the next nearest neighbor cou- 
pling are calculated by exact diagonalization in finite, two-dimensional Betts 
cells (lattices). A viable spin and charge coherent pairing, signaled by quan- 
tum critical points and the negative charge gap region, is found in 8- and 
10-site Betts lattices at small and moderate U regions consistent with our ex- 
act results in elementary bipartite geometries [Phys. Rev. B78, 075431 (2008)]. 
The contour isolines for continuous temperature driven-crossover between the 
Mott-Hubbard insulating and coherent pairing phases are demonstrated. The 
criteria for smooth and abrupt phase transitions are found for systematic en- 
hancement of coherent pairing by optimization of the next nearest neighbor 
coupling parameter. 

Keywords: high T c superconductivity, charge pairing, spin pairing, Betts 
lattice 
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1. Introduction 

A key element for understanding the physics in cuprates, pnictides, man- 
ganites and colossal magnetoresistant (CMR) nanomaterials is the experimen- 
tal observation of magnetic and density phase separation instabilities at the 
nanoscale The spatial inhomogeneities with short range electron corre- 

lations in the absence of long-range order play a crucial role in determining the 
mechanism of pairing instabilities in the high T c superconductors (HTSCs) and 
CMRs. Therefore, the local nature of the (pairing) interaction is sufficient for de- 
scribing electron pairing instabilities in the framework of the two-dimensional 
(2d) Hubbard model. The finite-cell-based lattice or small-size optimized clus- 
ters may be one of the few solid grounds available to handle this challenging 
issue A new, emerging guiding principle for the search of new materials 

can be identified as spatial inhomogeneities and density phase separation in- 
stabilities in the proximity to quantum critical points (QCPs). Strong quantum 
fluctuations can even dominate thermal fluctuations and affect the properties of 
a material well above absolute zero [4] (see Ref. [8] and also references therein). 
The phase sepa ration instabilities in small (bipartite) Hubbard clusters studied 
in Ref s. I9l ll0l,lllll display QCPs and interesting thermodynamics, which depend 
on the strength of the on-site Coulomb repulsion (U), the cluster topology and 
temperature. In contrast, spontaneous transitions in frustrated (non-bipartite) 
geometries can occur for all U by avoiding QCPs (level crossings) at finite IT 
and such transitions depend strongly on the sign of the hopping term f . 

The existence of the QCPs associated the phase separation instabilities may 
be a crucial ingredient of the superconducting transition and can provide im- 
portant clues for understanding the incipient microscopic mechanisms of charge 
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and spin pairing instabilities in the HTSCs and CMRs. Here our primary focus 
is the mechanism of electron pairing due to the density instabilities in small 
bipartite lattices in the ground state and finite temperatures, where there is no 
general agreement on the phase separation boundaries; related results are still 
controversial, especially at small U values fl2l Il3l, IT3l . Small systems suffer 
from finite-size and edge effects, so it is unclear whether the observed instabil- 
ities can survive in the thermodynamic limit. Thus, tests on reduced boundary 
effects are necessary to confirm the picture of instabilities in small cluster-based 
larger systems in the so called "optimized" Betts lattices 11511. 

The exact ground state properties of the infinite two-dimensional (2d) 
Heisenberg-like lattice have been extracted from systematic studies of finite 
Betts clusters 11611 . An infinite square lattice can be tiled by squares of L-site 
clusters with edge vectors I115ll which represent displacements of one vertex to 
the equivalent vertex in the neighboring tile with which it shares an edge (See 
Figure [l). Different block structures of "square symmetry" broken up into 
nearly decoupled clusters can be used as plaquettes to extrapolate the results 
to the infinite square lattice. Studies of finite clusters with periodic bound- 
aries and next nearest neighbor hopping can play an important part in testing 
the reliability of the results drawn from finite-size systems ll/Ulal . The phase 
separation boundaries have been calculated within the t-J model fl9ll20ll , al- 
though the validity of the strong-coupling expansion has been questioned 12111 . 
As far as the authors are aware, the exact studies of pairing thermodynamics 
have not been attempted even in small Betts-cluster-based lattices neither in the 
framework of single orbital Hubbard model with or without the next nearest 
neighbor coupling. 

In this paper, we address two specific issues: (i) check how exact calculations 
in finite Betts-cluster-based lattice are consistent with our exact results in 2x2 
and 2x4 clusters and, (ii) analyze the effect of lattice frustration on phase sep- 
aration instabilities in bipartite geometry driven by the next nearest neighbor 
coupling. A related question is whether the physical properties extracted from 
these Betts square lattices are likely to be applicable in the infinite lattice 11611 . In 
Sec. |2 the square geometry of the Betts lattice is explained and the correspond- 
ing Hamiltonian is described. In Sec. the exact results on charge and spin 
pairing in the Betts lattice based on canonical and grand canonical ensemble are 
analyzed and corresponding key physical properties are calculated. The finite 
temperature effects near QCPs at various U values are discussed in Sec. 13.1.41 
The enhancement of pairing driven by the next nearest neighbor coupling is 
studied in Sec 13.21 The conclusions in Sec. |4] provide a summary of what was 
learned from exact studies of electron pairing instabilities in the presence and 
absence of electron-hole symmetry. 



3 



©©©©©©©©© 
®©©©®©©®© 
© © © © © ® © © © 
©©® ©©®® ©© 
© © © © © © © © © 
© © © © ® <2K© © © 
® © © © (2(© © © © 
© © ® © © © © 
© © © © © ©%> © © 

Figure 1: The 8-site finite unit cell (plaquette) for the square lattice. When repeated periodically, it 
can fill the entire (infinite) space. The cells have li and I2 edge vectors, (2,2) and (2,-2), as defined 
inRef. (H. 



4 



2. Betts square lattice 

As shown in Figured] the entire 2d lattice space can be tiled by periodically 
repeated small Betts squares, which are connected through edge vectors li and 
I2. Such periodic boundary conditions reduce the edge effects while taking 
full advantage of the space group symmetries of the isotropic two dimensional 
(square) lattice. Notice that, in the 8-site Betts lattice, each odd site is surrounded 
by all the even sites (as nearest neighbors) and vice versa. This structure 
illustrates the bipartite character of Betts (square) cells. For comparison, we 
study the effect of frustration on electron pairing in the absence of electron-hole 
symmetry by introducing the next nearest neighbor coupling t nnn with periodic 
boundary conditions. The model describes the electron doping for t mm > and 
the hole doping regime for t nnn < 0. 

2.1. Model 

We consider the extended Hubbard model 



where summation over i and / in Eq. (T) goes through all lattice sites L with 
coupling integral tu 



and LI > is the on-site Coulomb interaction. The energies are measured with 
respect to t nn > 0, which is set to 1 everywhere, unless otherwise stated. 




(1) 




(2) 
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Figure 2: Schematic drawing of hole (electron) redistribution within Hubbard nanoclusters with 
one hole of half-filling in the grand canonical ensemble of 8-site clusters at low temperatures. The 
state on the left is an ordinary state at N = 7 per unit cell. Notice that the spontaneous fluctuations 
in particle number, near the average (N) = 7, are energetically favorable and make electron redis- 
tribution across the ensemble of clusters possible even without direct contact (hopping) between 
them. When A c > 0, the hole localization on separate clusters corresponds to a homogeneous 
stable Mott-Hubbard-like insulating tf state with electron-hole pairing. When A c < 0, the hole pair 
binding in d 6 and half-filled d 8 inhomogeneous cluster configurations on the right are energetically 
preferred. 
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3. Results 



3.2. Phase separation instabilities 

We consider L=8 and 10 site Betts clusters where the dimensions of the 
Hilbert spaces are 4 s and 4 10 respectively. An exact diagonalization along with 
the Lanczos algorithm [22] is employed to evaluate the relevant, low lying 
eigenstates at all filling (doping) values. These eigenstates are used in Sec. |3.1.H 
and l3.1.3| to extract the charge and spin pairing gaps and corresponding pairing 
instabilities of the model. 



3.2.2. Charge pairing instability 

The charge pairing instability in the negative charge gap region, as defined 
below, is a precursor to a spin pairing instability at rather low temperatures 
(see Sec. 13. 1.3| |. In our aforementioned publications, we have discussed this 
issue in selected cluster geometries such as the 2x2 square and 2x4 ladder. 
We identify charge pairing behavior at zero and finite temperatures in a L-site 
cluster by defining a charge gap A C (N, T) at a given U in a particular doping 
region: 

A C (N, T) = E(N + 1, T) + E(N - 1, T) - 2E(N, T) (3) 

where E(N, T) is the lowest canonical many body energy for an N-electron 
state at a fixed temperature T. This charge excitation gap determines the 
stability of an N— electron state compared to an equal admixture of (N + 1)— 
and (N - 1)— electron states. As shown in Figure [2[ the electron number in 
clusters placed in a thermal reservoir can fluctuate near an average number (N). 
This schematic picture for (N) = 7 illustrates the incipient mechanism of hole 
pairing at the 8-site cluster level, where properties such as density can strongly 
fluctuate. (These fluctuations closely resemble interconfiguration fluctuations 
in mixed valence compounds j23ll .) Eq. © describes an interconfiguration 
energy gap for electron fluctuations between different many body states and 
cluster configurations (d) that differ in electron number N [24], 

d N + d N = d N+l +d N -\ (4) 

Depending on the strength of the on-site electron-electron repulsion U, this 
charge excitation gap can be positive A c > or negative A e < 0. Physically, 
A c > manifests a stable d N (cluster) configuration, while A c < describes a first 
order phase separation instability with spontaneous generation of d N ~ x and d N+1 
electron (cluster) configurations. Thus, a negative charge gap (A c < 0) implies 
an effective attraction or pairing between holes accompanied by a (spontaneous) 
transition into an inhomogeneous phase, while a positive A c shows a preference 
for holes to be localized in separate clusters which are intimately related to 
insulating behavior and the Mott-Hubbard-like electron-hole pairing [24]. For 
a given chemical potential and temperature in the grand canonical ensemble, 
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Figure 3: The canonical charge A c and spin A s excitation gaps in the ensemble of the 8-site Betts 
cluster for N - 7 as a function of U at f„„„ = and infinitesimal T —> 0. The grand canonical spin 
gap A* > provides the stability for spin singlet pairing at < U < U c , The nodes A c ' s ' h (U) = 
of the charge and spin gaps define a quantum critical points for various electron instabilities. 
This phase diagram displays the main phases found earlier for elementary square geometry in 
Ref. [21. Phase I is a phase with a negative charge gap and a positive spin gap of equal amplitude 
(A c = -A ,! ) at U < U c which describes coherent pairing of electrons (holes) in the spin singlet 
state. Phase II is a spin liquid phase with gapless, low spin- 1 excitations separated from the higher 
spin- 1 by canonical gap A s > 0. In Phase III, the negative spin gap manifests the onset of low 
spin- 1 (unsaturated) ferromagnetism. (U u identifies the crossover point for the spin-z 3/2 to 5/2 
transition.) There is also a consequent transition into a fully saturated spin-| ferromagnetism at 
larger a U value, U s (not shown). 
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Figure 4: A plot of the charge gap A c in the 10-site Betts cell at N = 9 as a function of U at t nnn = 
and T = 0. In contrast to the 8-site lattice, optimized pairing is shifted to larger U values and a 
positive charge gap appears at very small U values. This is due to the reduced symmetry of 10-site 
geometry (see the inset) compared with the 8-site lattice (see Sec [3X2) . 



electrons are allowed to be redistributed among the clusters to optimize their 
free energy. Thus, on the left of Figure a homogeneous Mott-Hubbard 
state with electron-hole pairing is preferred in the positive charge gap (A c > 0) 
region, while in the negative gap region (A c < 0) on the right, an inhomogeneous 
phase with hole-rich (N = 6) and hole-poor (N = 8 half-filled) clusters becomes 
energetically favorable. 

The charge gap A c at T = is plotted as a function of U in Figure [3] for 
N = 7. The nodes for the sign change in the excitation gap provide strong 
evidence for the existence of QCPs for a level crossing instability associated 
with a first order phase separation transition [9]. For example, the charge gap 
in Figure [3]vanishes (A c (ii c ) = 0) at IT = and also at a particular U value, i.e., 
U c = 8.52. These QCPs describe an abrupt (discontinuous) change in the ground 
state of a many-body system due to strong quantum fluctuations. Compared 
to the elementary square geometries, the negative charge gap (A c < 0) for 
the 8-site clusters with electron number N = 7 has a larger (antinode gap) 
maximum binding energy and shows a larger region < U < U c of electron 
instability. The Mott-Hubbard insulating ground state for electron-hole pairing 
with A c > is stable at all U > U c (shown in Fig. [3) while for < U < U c , the 
ground state consists of paired electrons. This instability in redistribution of 
electrons, caused by electron hopping between clusters, signals the formation 
of inhomogeneous, hole-poor d N ~ l and hole-rich d N+1 regions or indicates a 
tendency toward phase separation |24Tj . 



3.1.2. Square symmetry 
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Broken symmetry of the two-dimensional planar geometry in various finite- 
size Hubbard clusters plays a crucial role in determining the pairing symmetry 
and superconducting properties. Our results in various square (bipartite) ge- 
ometries show that pairing properties strongly depend on the cluster symmetry. 
The gap behavior in the 8-site Betts cells is different from the 2x4 ladder jltjllTll . 
For comparison, in Figure |4j we show the exact charge gap with reduced square 
symmetry in the 10-site Betts lattice as a function of IT. The shallow charge gap 
at small U is transformed to a significantly larger one with maximum (antin- 
ode) binding gap energy at moderate U. This figure also shows an oscillating 
gap due to the reduced square symmetry which resembles that of the 2x4 
clusters yjj. Both the 2 X 4 ladder and the 10-site Betts lattice (with edge 
vectors (1,3) and (3,-1)) have lower symmetries than the 8-site Betts cell (with 
edge vectors (2,2) and (2,-2)) and we believe that the oscillations in the charge 
gap (as a function of U) are related to this. Although a positive charge gap 
appears at very small IT values, in overall, the evolution of the gap is quite 
similar to our results obtained in other clusters, except at very low IT. Thus, 
Betts 8- and 10-site clusters provide strong support for electron instabilities and 
nanoscale inhomogeneities found in generic 2x2 and 2x4 clusters in Gum 
and reproduced later [2g]. 

3.2.3. Spin (pseudo)gaps 

Spin (pseudo)gaps classified according to their total spin S give essential 
insights into the behavior of the system on response to applied magnetic field 
h in the vicinity of one hole off half filling for negative (IT < U c ) or positive 
(IT > U c ) charge gaps. This is simply related to a ground state having total spin 
either zero or nonzero respectively. 

Using exact analytical expressions for the grand canonical potential we analyze 

the variation of the spin S z (p) and susceptibility Xh — —mr as a function of 
magnetic field h near the critical chemical potential, p c , close to given average 
electron number, (N) = L — 1. In the negative charge gap region, near u c , the 
ground state has zero total spin, S = 0. In order to break spin pairs and generate 
spin excitations, a finite magnetic field need to be applied. The threshold value 
of this field will be identified as a spin pseudogap (A ,! ), as described in our 
previous publications |@,|l3]. This spin pseudogap at < U < U c provides 
(rigidity) stability (A'' > 0) to paired (singlet) spins in the ground state. At rather 
low temperatures, near T — > 0, the stable spin pseudogap is equal in magnitude 
to the charge pairing gap, A p = -A c . This is an important characteristic of phase 
coherence in the ground state associated with simultaneous coherent pairing 
of independent charge and spin (Fermion) entities and their corresponding 
full Bose-Einstein condensation (real-space pairing) with a single energy gap. 
Such behavior is similar to the coherent pairing with unique quasiparticle gap 
in conventional BCS theory [9, 24]. It is important to note that the coherent 
pair formation through separate condensation of bound holes and coupled 
opposite spins (both bosons) at different critical temperatures T c p and T s p in 
real space here is somewhat different from the Cooper pairs (composite bosons) 
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condensation in momentum space with the unique critical temperature. The 
magnitudes of the spin pseudogap and charge gap also differ with increasing 
temperature here, unlike in the BCS case. It is consistent with ARPES data that 
excitations in high T c superconductors (as in the exact solution) are not well 
defined quasiparticles as it is in the BCS theory [8]. 

In the positive charge gap U region(s), the behavior of states near \i c is 
different from the above due to the presence of a spin-charge separation effect 1 9] 
in the ground state that carries a nonzero total spin, S + 0. Clearly, this ground 
state can be excited with an infinitesimal magnetic field, hence spin excitations 
here are gapless. In addition to the pseudospin gap A/,, one can also focus on 
the total spin S of the ground state which has a spin degeneracy of (2S+1). 
It turns out that, depending on the value of U due to level crossings, certain 
competing states with high or low (total) spin can become the ground state. 
Hence, we define a canonical spin gap A s as the energy difference between the 
ground state with spin S and the lowest excited state with spin S having the 
same N and U as 

_ f E(N, S' r T) - E(N, S, T) when S' > S, 
( ' ' \ E(N, S, T) - E(N, S', T) when S' < S, { ' 

where E(N, S, T) is the canonical energy of an N electron state with total spin 
S at rather low temperature T. If the excited state has a higher (lower) total 
spin, the spin gap is positive (negative). This gap (along with its sign) identifies 
whether a high or low spin state turns out to be the lowest average energy 
state (or ground state at T — > 0). In addition, the ground state here can be 
described as a spin liquid due to the spin degeneracy (2S+1), which can be 
lifted by infinitesimal h — > 0. One other noteworthy point is that, in the 8-site 
Betts lattice, there are several U values (QCPs) at which such level crossings 
occur (from low-spin to high-spin or vice versa) near \i c . Those are denoted by 
U c (spin to 3/2), U u (from spin 3/2 to 5/2) (see Fig. |3]l and U s (from spin 5/2 to 
7/2). 

3.2.4. Finite temperature effects 

The QCPs introduced in Sees. 13.1. H and I3.1.3l turn out to be useful for the 
analysis of spontaneous instabilities at nonzero temperatures due to the inter- 
play of quantum and thermal fluctuations in proximity to the above discussed 
QCPs. These QCPs as well as the doping dependencies on the chemical po- 
tential at nonzero crossover temperatures (not shown here) are qualitatively 
similar to those obtained earlier [24]. A contour line (isoline) for the charge 
gap as a function of T and If, along which the gap attains a constant value c 
(A C (U, T) = c) defines the contour map in Figure[5] The slice cut of the vanishing 
gap (c = 0) at finite temperatures (A C (U, T) = 0) defines the boundary between 
positive and negative charge gap phases, which is marked by the solid line in 
Figure [5] This crossover is a smooth second order transition for the onset of 
pairing instabilities at finite temperatures. For a given T, the magnitude and 
area of the negative gap region are an indication of the pairing strength and the 
region of the phase separation instability; this magnitude and area shrink and 
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Figure 5: A contour plot of charge gap in the 8-site (periodic) Betts cluster at N = 7 as a function 
of Coulomb repulsion U and temperature T. As T increases, the region of negative charge gaps 
becomes smaller and the crossover point U C (T) shifts to lower U values. The negative gap disap- 
pears completely at about T = 0.16. The solid boundary in the figure denotes the contour (U) 
at which the charge gap vanishes (A c = 0). Phase I with a shaded pattern is the electron-electron 
pairing phase. Phase II without a shaded pattern is the electron-hole pairing phase. 

disappear at relatively low temperature, indicating that the pairing and phase 
separation instability is a low temperature phenomenon. The solid boundary 
line evolves smoothly with temperature and approaches the QCPs, U = and 
U c = 8.52, as temperature goes to zero, T — > 0. Exactly at T = 0, this picture 
shows the region < U < U c for coherent pairing. The negative gap region 
defines first order phase separation instabilities, while the positive gap is a sig- 
nature for smooth second order transitions. Notice that the electron pairing gap 
below T C P {U) with A C (T) < is transformed into a Mott-Hubbard gap A C (T) > 
when the temperature increases above T C P (U). This picture describes the tran- 
sition from incoherent charge pairing of preformed electron (hole) pairs with 
decoupled spin into a Mott-Hubbard insulator. Thus, the continuous transition 
at the classical critical temperatures here from one, ordered, low-temperature 
phase into another (quantum) ordered high temperature phase is somewhat 
different from the conventional phase boundary described by an order param- 
eter, which is non-zero in the ordered phase and zero in the disordered high 
temperature phase. Note, that the high temperature insulating phase in Fig- 
ure [5] at T C F (U) in < U < U c region originates from the low temperature 
(quantum) Mott-Hubbard phase at U > U c . Thus, we conclude that the QCP 
here is unique, since it describes a boundary for a smooth transition between 
two quantum ordered phases with different symmetries. The effect of quantum 
criticality in Figure|5]can be felt at low temperatures without even ever reaching 
the ground state (T = 0). 



12 



3.2. Next nearest neighbor coupling 

In many real materials (such as the cuprates), the contribution of hop- 
ping among next nearest neighbor (ram) atoms can be important due to many 
reasons, such as its ability to separate hole and electron doped properties by 
breaking particle-hole symmetry 12611 . Experimental results strongly support 
the idea that holes in CuC>2 planes are mainly sited in the d x i_ y 2 Cu orbital in re- 
lated undoped insulating materials, and additional holes are transferred mainly 
to the oxygen p orbitals upon doping. The next-nearest coupling term has al- 
ready been invoked to explain the shapes of the Fermi surfaces obtained from 
recent band-structure calculations for these new superconducting materials. 
The two orbital Hubbard model with p-bonding orbitals, after elimination of 
the Cu sites, can be reduced to a single-orbital Hubbard model with nearest and 
next-nearest-neighbor couplings 10 by allowing holes to move within a given 
oxygen sublattice. The pairing can be changed significantly by the Coulomb 
repulsion, when the second-neighbor interactions are included. Therefore, a 
study in the Hubbard model 10 with t nnn + might give a more realistic phys- 
ical picture related to real materials. In this section, we discuss the effects due 
to a nonzero t nnn . Because t mm in real material is usually much smaller than t 
parameter we only show results for -0.3 < t nnn < 0.3. Below we will find con- 
ditions under which broken C4 symmetry in frustrated Betts cells with t nnn + 
can be harmful or favorable for pair binding. 

At zero temperature, charge gaps at one hole off half- filling (contour plots or 
isolines) for different t nnn and LI values are shown in Figure |6] The charge gap 
behavior in the region —0.3 < t nnn < 0.3 is similar to that found in section 13.1.11 
A negative charge gap, representing charge pairing, can be found at relatively 
small IT values. The bold line in Figure |6]shows where the charge gap becomes 
zero, which defines a boundary between charge pair phase and electron-hole 
pairing phase. As t nnn increases from negative (opposite sign as that of t) to 
positive (same sign as that of f), the boundary shifts to larger U values and the 
magnitude of the gap also becomes larger. This can be explained by the effects 
of t mm on the ground state with holes. The hopping between next nearest 
neighbors are energetically unfavorable for the half-filled antiferromagnetic 
state since the next nearest neighbors have the same spin direction. When 
holes are created in the antiferromagnetic background, hopping between next 
nearest neighbors will become more active and will have a measurable effect on 
the ground state. Depending on the sign of t nnn , it will make the ground state 
with holes higher (t nnn < 0) or lower (t nnn > 0) in energy. Therefore a positive 
t nnn makes hole-hole pairing (charge pairing) ground state more energetically 
favorable and enhances the negative charge gap. In spite of this introduced 
frustration in the square system, pairing can be enhanced by appropriate sign 
of the next-nearest-neighbor coupling (£„„„ > 0). 

At fixed U, tracing charge gaps for different t nnn yields Figure [7] In the 
tnnn < region, charge gaps are linear with respect to t mm , while they show 
nonlinear behavior for t mm > 0. The charge gap for T = is consistent with the 
conclusions drawn from ground state calculations in t — t' — } model at t' < 
and f > |2^|. This linear behavior for negative t nnn suggests that the variation 
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Figure 6: A contour plot of charge gap A c for different U and t nnn at T - 0. As t nnn increases the 
crossover point shifts to a larger U values and the maximal charge gap width increases as t nnn > 
but remains the same as t nnn < 0. The negative gap region is shaded. 

of the charge gap for small nonzero \t nnn \ can be obtained from the case for zero 
tnnn with a linear extrapolation, which makes further research on next nearest 
hopping easier. Further calculations show that the linear behavior is lost when 
tnnn is positive and large; since t nnn has already exceeded the range where t nnn is 
reasonably small (as the hopping between next nearest neighbors), we will not 
discuss it in this article. Figure [7] also shows that the linear slope approaches 
zero as the magnitude of the gap approaches its maximum (antinode) value 
(around LI = 4), which indicates that, although negative t nnn drives charge 
pairing crossover point to smaller U values, it does not affect the maximum 
gap value significantly. 

4. Summary 

We have discussed charge and spin pairing in the 8- and 10-site Betts 
clusters subject to periodic boundary conditions. This work provides ample 
evidence that such pairing instabilities do exist and are robust at the 8- and 10- 
site cluster sizes in the ground state and at finite temperatures. An important 
question is whether the obtained electron instabilities will continue to exist 
in the two-dimensional lattice as the cluster size increases, especially, in the 
thermodynamic limit. As exact calculations cannot go up to clusters large 
enough to eliminate size effects, the periodic Betts cells are considered to be 
the best optimal structures that can minimize and reduce edge effects. The 
key intrinsic properties of an infinite square lattice can be extracted from exact 
calculations in finite Betts lattices B17I1 . For example, the extrapolated ground 
state energy from 8- and 10-site Betts cells per site at half-filling at large U 
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Figure 7: Charge gap at fixed U value for different t nnn . A linear relation is found in the negative 
from region. When t„„„ > 0, the variation of the gap is not linear. This picture suggests that clusters 
with different signs for the next nearest neighbor hopping have different effects on the spectral 
properties and electron charge pairing. 

is quite close to existing analytical result in thermodynamical limit, —1.15 X 
l28ll . The basic scenario of electron pairing due to the phase separation 
instabilities, reproduced in small Betts lattice near half-filling, can be valid 
also in larger size Betts lattices. Our finite size studies are also quite relevant in 
view of the observed magnetic and density phase separation instabilities |H|] and 
electronic inhomogeneities at the nanoscale in cuprates, pnictides, manganites 
and colossal magnetoresistance (CMR) nanomaterials [24]. 

The exact diagonalization of the Hubbard model on Betts cells (near one 
hole off half -filling) can be used to extract QCPs and crossover temperatures, 
which in turn provide important insight into the mechanisms of pairing and 
phase separation instabilities I9l 12911 . The exact diagonalization provides strict 
criteria for phase separation instabilities and boundaries between smooth and 
first order phase transitions [23]. The introduction of next nearest neighbor 
hopping shifts the quantum crossover points but cannot eliminate the condi- 
tions necessary for electron and spin pairing. Moreover, we find that by an 
appropriate choice of the sign of t nnn > 0, electron charge and opposite spin 
pairings in coherent phase can be strongly enhanced. The sign change of t nn „ 
term have different effect on the spin pairing. In contrast, we find a stabiliza- 
tion of the Nagaoka ferromagnetism with the negative spin pairing gap of the 
optimal value at large LI values in case of negative sign of t nn „ < 0. 
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